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1 Introduction 



The stochastic heat equation 
d v d 2 



u(t,x) = p(u(t,x))W(t,x), x G R, t G 



Ot 20x 2 J" K " , ~ J n " v "'"" " "^~+> (l.i) 

u(0,-) = /i(-) , 

where is space-time white noise, p(u) is globally Lipschitz, is the initial data, and 
E* = ]0, +00 [, has been intensively studied during last two decades by many authors [1, 2, 
3, 4, 5, 7, 8, 9, 10, 11, 13, 12, 14, 15] . In particular, the special case p{u) = \u is called the 
parabolic Anderson model [2]. Our work focuses on this equation with general deterministic 
initial data, and we study how the initial data affects the solution. 
The one-dimensional heat kernel function is 



G v (t,x):= expt-^p-}, (t,x)eR* + xR. (1.2) 



12 

y/2nut ""^ I 2z4 

For the existence of random field solutions to (1.1), the case where the initial data p is a 
bounded and measurable function is covered by the classical theory of Walsh [15]. When p 
is a positive Borel measure on R such that 

sup sup Vr (p, * G v (t, o)) (x) < 00, for all T > 0, (1.3) 

te[o,T] xm 

where * denotes convolution in the spatial variable, Bertini and Cancrini [1] gave an ad- 
hoc definition for the Anderson model via a smoothing of the space-time white noise and a 
Feynman-Kac type formula. Their analysis depends heavily on properties of the local times of 
Brownian bridges. Recently, Conus and Khoshnevisan [5] constructed a weak solution defined 
through certain norms on random fields. The initial data has to verify certain technical 
conditions, which include the Dirac delta function in some of their cases. In particular, 
the solution is defined for almost all (t,x), but not at specific (t,x). More recently, Conus, 
Joseph, Khoshnevisan and Shiu [3] also studied random field solutions. In particular, they 
require the initial data to be a finite measure of compact support. We improve the existence 
result by working under a much weaker condition on initial data, namely, p can be any 
signed Borel measure over R such that 

(\p\ *G v (t, -))(x) < +00 , for allt > and x G R , (1.4) 

where, from the Jordan decomposition, p = p + — p_ where p± are two non-negative Borel 
measures and \p\ := p + + //_. For instance, if p(dx) = f(x)dx, then f(x) = exp(a|x| p ), 
a > 0, p G ]0,2[, (i.e., exponential growth at ±00), will satisfy this condition. Proposition 
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2.9 below shows that initial data cannot be extended beyond measures to other Schwartz 
distributions, even with compact support. 

Moreover, we obtain estimates for the moments E(|tt(t, x)\ p ) with both t and x fixed for 
all even integers p > 2. In particular, for the parabolic Anderson model, we give an explicit 
formula for the second moment of the solution. When the initial data is either Lebesgue 
measure or the Dirac delta function, we give explicit formulae for the two-point correlation 
functions (see (2.21) and (2.24) below), which can be compared to the integral form in Bertini 
and Cancrini's paper [1, Corollaries 2.4 and 2.5]. 

Our proof of existence is based on the standard Picard iteration. The main difference 
from the conventional situation is that instead of applying Gronwall's lemma to bound the 
second moment from above, we show that the sequence of the second moments in the Picard 
iteration converges to an explicit formula (in the case of the parabolic Anderson model). 

After establishing the existence of random field solutions, we study whether the solution 
exhibits intermittency properties. More precisely, define the upper and lower Lyapunov 
exponents for constant initial data (Lebesgue measure) as follows 



V = i imsup !2^M^ffi, v=liminf !2iiM^M (L5) 

t^+oo t P t^+oo t 

Following Bertini and Cancrini [1], we say that the solution is intermittent if A n := A n = A n 
and the strict inequalities 

*,<£<...<£<••. (1.6) 
2 n 

are satisfied. Carmona and Molchanov gave the following definition [2, Definition III. 1.1, on 
p. 55]: 

Definition 1.1 (Intermittency). Let p be the smallest integer for which X p > 0. When 
p < oo, we say that the solution u(t, x) shows (asymptotic) intermittency of order p and full 
intermittency when p — 2. 

They showed that full intermittency implies the intermittency defined by (1.6) (see [2, 
III. 1.2, on p. 55]). This mathematical definition of intermittency is related to the property 
that the solutions develop high peaks on some small "islands". The parabolic Anderson 
model has been well studied: see [2] for a discrete approximation and [1] and [9] for the 
continuous version. Further discussion can be found in [16]. 

When the initial data are not homogeneous, in particular, when they have certain ex- 
ponential decrease at infinity, Conus and Khoshnevisan [4] defined the following lower and 
upper exponential growth indices: 

A(n) :=sup < a > : limsup- sup logE(|u(t, a;)| n ) > > , (1.7) 

I <->oo t \ x \> a t \ 
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X(n) :=inf <( a > : limsup- sup logE (\u(t, x)\ n ) < } , (11 

t^tco t \ x \> at 

and proved that if the initial data p is a non-negative, lower semicontinuous function with 
compact support of positive measure, then for the Anderson model (p(u) = Xu), 

f < A(2) < A(2) < f . 

We improve this result by showing that A(2) = A(2) = A 2 /2, and extend this to more general 
measure-valued initial data. This is possible mainly thanks to our explicit formula for the 
second moment. 

We now discuss the regularity of the random field solution. Denote by C^p 2 (D) the 
set of random fields whose trajectories are almost surely /3i-H61der continuous in time and 
/^-Holder continuous in space on the domain D C R + x R, and let 

C^-ip) := |J |J C aua2 (D). 

C*lS ]0,y8l[«26 ]0,/3 2 [ 

In Walsh's notes [15, Corollary 3.4, p. 318], a slightly different equation was studied and 
the Holder exponents given (for both space and time) are 1/4 — e. Bertini and Cancrini [1] 
stated in their paper that the random field solution for the parabolic Anderson model with 
initial data satisfying (1.3) belongs to Ci_ i_(M^_ x R). In [11, 14], the authors showed that 
if the initial data is a continuous function with certain exponentially growing tails, then 



u e Ci_ i_(R+ x R), a.s. (1.9) 

4 '2 

Sanz-Sole and Sarra [13] considered the stochastic heat equation over M. d with spatially 
homogeneous colored noise which is white in time. Let p, be the spectral measure satisfying 



/ T\ i \c\2\v < + 00 ' for some V e ]°> ![■ 

JR d I 1 + Is I J 



They proved that if the initial data is a bounded p-H61der continuous function for some 
p G ]0, 1[, then the solution is in 

u(t,x) G Ci (pA{1 _^ )) _ !pA{1 _ 7?) _ (M* + x R) . 

For the case of space-time white noise, the spectral measure jx is Lebesgue measure and 
hence r] can be 1/2 — e for any e > 0. Their result ([12, Theorem 4.3]) reduces to 

tt(t,s)6C (iAf)»-^ xR ) • 
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More recently, Conus et al proved in their paper [3, Lemma 9.3] that the random field 
solution is Holder continuous in x with exponent 1/2 — e (for initial data that is a finite 
measure). They did not give the regularity estimate over the time variable. In their pa- 
pers [7, 8] , Dalang, Khoshnevisan and Nualart considered a system of heat equations with 
vanishing initial conditions subject to space-time white noise, and proved that the solution 
is jointly Holder continuous with exponents 1/4— in time and 1/2— in space. We extend 
the Ci i_ (R* x R) -Holder continuity result to measure-valued initial data satisfying (1.4). 
We show that the result in (1.9) should exclude the time line t = unless the initial data 
is l/2-H61der continuous. 

The difficulties for the proof of the Holder continuity of the random field solution lie in 
the fact that for the initial data satisfying (1.4), the p-th moment E [\u(t, x)\ p ] is neither 
bounded for iGl, nor for t G [0, T]. Standard techniques, which isolate the effects of initial 
data by the L p (£7)-boundedness of the solution, fail in our case. Instead, the initial data play 
an active role in our proof. Note that Fourier transforms are not applicable here because \i 
need not be a tempered measure. 



2 Main Results 

Denote the solution to the homogeneous equation 

d v d 2 



u(t,x) = 0, x G R, t G 



dt 2&V ' ' +' (2.1) 

u(0, •) = fi(-) , 

by 

J (t, x) := fa * G u (t, •)) (x) = / G v (t, x - y)n(dy) , (t, x) G MI x R . 



Note that Jo(t,x) is well-defined by the hypothesis (1.4). We formally rewrite the stochastic 
partial differential equation (1.1) in the integral form (mild form): 

u(t,x) = J (t,x) + I(t,x) (2.2) 

where 



I(t, x) := G v {t-s,x-y)p (u(s, y)) W(ds } dy) . (2.3) 

J J [0,t]xR 

By convention, 1(0, x) = 0. The above stochastic integral is defined in the sense of Walsh 
[15, 6]. 
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2.1 Notations and Conventions 

Assume that the function p : R i— > R is globally Lipschitz continuous with Lipschitz constant 
Lip p > 0. We need some growth conditions on p: Assume that 

\p{x)\ 2 < L\ +|x| 2 ) , for all x e R , (2.4) 

for some constants L p > and q > 0. Note that \/2Lip p < L p , and the inequality may be 
strict. In order to bound the second moment from below, we will sometimes assume that for 
some constants l p > and £ > 0, 

\p(x)\ 2 >l 2 p (q 2 +\x\ 2 ) , forallxGM. (2.5) 

We shall also give special attention to the linear case (the parabolic Anderson model) p(u) = 
Xu with A ^ 0, which is a special case of the following quasi-linear growth condition: for 
some constant q > 0, 

\p(x)\ 2 = X 2 (<? 2 + |x| 2 ) , for all x e R . (2.6) 



We use the convention that G u (t, •) = if t < 0. Hence, the integral region in the 
stochastic integral in (2.3) can be written as R + x R. 
Define a kernel function 

K { tM:=G s bx)(^ + £e£*(*J±) ] , (2.7) 

for all (t,x) G W* + x R, where <&(x) is the probability distribution function of the standard 
normal distribution: 

r x e -v 2 / 2 
^ := / ~^^ dy ■ 

J-oo V^TT 

We also use the error function erf(x) := j^e~ y2 dy and its complement erfc(x) := 1 — 
erf(x). Clearly, 



$(x) = i (l + erf (x/v^)) , 
erf (x) = 2$ (\/2 a;) - 1, erfc(x) = 2 (l - $ f\/2 ; 
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We use * to denote the simultaneous convolution in both space and time variables. Define 
another function 



n(t; is, A) := (1 * K) (t, x) = 2e^$ A ; 



t 




2v 



1 . 



(2.8) 



Clearly, /C (t, x; v, A) can be written as 

K (t,x; v, A) = G u/2 (t,x) 
We use the following conventions: 



( \jAixvt ' 4z/ 



(2.9) 
(2.10) 
(2.11) 

p ) , forallp>2, (2.12) 

where z p is the universal constant in the Burkholder-Davis-Gundy inequality (in particular, 
z 2 = 1) and (Xp £ IS db constant defined as 



K,(t, x) 


:= K.(t,x ; 


z/,A) , 


JC(t, x) 


:= K (t, x ; 


^,£p) > 


tC(t, x) 


:= K.(t,x ; 


M P ) , 


JC p (t,x) 


:= K (t, x ; 





'2(p-V/p if^O, p>2, 
V2 if ^ = 0, p>2, 
1 ifp = 2. 



(2.13) 



We only need to keep in mind that a p ^ < 2. Note that the kernel function JC p (t, x) implicitly 
depends on q through a p ^ which will be clear from the context. If p = 2, then /C 2 (t,2;) = 
JC(t, x). 

Similarly H(t), H(t) and H p (t) denote the kernel functions with A in H(t) replaced by 
L p , l p and a p ^z p L p , respectively. Again H p (t) depends on q implicitly which will be clear 
from the context. 

Let us set up the filtered probability space. Let 



| W t (A) : AeB b (R),t >o} 



be a space-time white noise defined on a probability space (£l,J-, P), where B b (R) is the 
collection of Borel measurable sets with finite Lebesgue measure. Let (Tt,t > 0) be the 
standard filtration generated by this space-time white noise. More precisely, let 



jro ;= a (W S {A) : < s < t, A e B b (R)) V A/", t>0 
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be the natural filtration augmented by all P-null sets in T . Define JF t := = A s>i J r ° 
for any t > 0. In the following, we fix this filtered probability space {O, J 7 , {J-j : t > 0}, P}. 
We use 1 1-|| to denote the L p (f2)-norm. Denote \p~\2 '■= 2|~p/2], which is the smallest even 
integer greater than or equal to p. 

Let Ai (R) be the set of locally finite Borel measures over R. Define 

Mg(M) := |/i G M(R) : ^ e^ |a:| |//|(dx) < +00 J , /3 > 0, (2.14) 

where = p + + /U_ is the Jordan decomposition of a measure into two non- negative 
measures. We use subscript "+" to denote the subset of non-negative measures. For example, 
Ai + (M.) is the set of non-negative Borel measures over MandA4^ + (R) = A4^(R)nA4+(R). 

A random field Y = (Y(t, x) : (t, x) G R* + x K) is said to be L P (Q) -continuous, p > 2, if 
for all (t, x) G K* x K, 

lim =0. 

2.2 Existence, Uniqueness and Moments 

We first give the definition of the random field solution as follows: 

Definition 2.1. A solution u = (u(t,x) : (t, x) Gl^xl) to (1.1) (or (2.2)) is called a 
random field solution if 

(1) it is adapted, i.e., for all (t,x) G R* + x R, u(t,x) is J-^-measurable; 

(2) u is jointly measurable with respect to i3 (R* + x R) x J 7 ; 

(3) (G^ * I \p(u) 1 12) (t,x) < +00 for all (t,x) G R* + x E, and the function (t,x) ^ J(t,x) 
mapping from Kl x R into L 2 (Q) is continuous; 

(4) m satisfies (1.1) (or (2.2)) almost surely, for all (t,x) G R* + x R. 
The first main result is stated as follows. 

Theorem 2.2 (Existence, uniqueness, moments). Suppose that 

(i) the initial data p, is a signed Borel measure such that (1.4) holds; 

(ii) the function p is Lipschitz continuous such that the linear growth condition (2.4) holds. 

Then the stochastic integral equation (2.2) has a random field solution u = {u(t,x) : t > 
0,x G R} (note that t > 0) in the sense of Definition 2.1. This solution has the following 
properties: 
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(1) u is unique (in the sense of versions); 

(2) (t,x) i — y u(t,x) is L P (Q)- continuous for all integers p > 2; 

(3) For all even integers p > 2, the p-th moment of the solution u(t,x) satisfies the upper 
bound 

f J 2 (t,x) + (J 2 */C) (t,x) %{t), ifp = 2, 

\\u(t,x)\\l<\ (2.15) 
\2J 2 (t,x) + (2J 2 *£ P ) (t,x) + ^ 2 U p (t), ifp > 2, 

for all t > 0, x G IR, and the two-point correlation satisfies the upper bound 

E[u{t,x)u{t,y)} 

< J (t, x)J (t, y) + J 2 / ds f(s, z)G v {t -s,x- z)G v {t - s,y - z)dz 



+ ^\x-y\($ (^j=4 )-l)+2L 2 p <?t G 2v {t, x-y), (2.16) 



<2vt 

for all t > 0, x,y G IR, where f(s, z) denotes the right hand side of (2.15) forp = 2; 
(4) If p satisfies (2.5), then the second moment satisfies the lower bound 

\\u{t,x)\\l > J 2 (t,x) + (J 2 */C) (t,x) +S 2 H(t) (2.17) 
for all t > 0, x G IR, and the two-point correlation satisfies the lower bound 

E[u(t,x)u(t,y)} 

> J (t, x) J (t, y) + l 2 p ds f{s, z)G v it - s, x - z)G v {t - s,y - z)dz 



+ l i^\x - y\ ( $ (%=4 ) 1 ) + 2 lW t G 2u (t, x-y), (2.18) 



/2i/t 

/or a// 1 > 0, G IR, where f(s,z) denotes the right hand side of (2.17); 

(5) In particular, for the quasi-linear case \p{u)\ 2 = A 2 (<; 2 +u 2 ), the second moment has the 
explicit expression 

\\u(t,x)\\ 2 2 = Jl(t,x) + (Jo**:) (t,x) + C 2 Hit) , (2.19) 
/or all t > 0, x G IR, and £/ie two-point correlation is given by 
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E [u(t, x)u (t, y)] 



o 



J (t, x)J (t, y) + A ds f(s,z)G v (t — s,x — z)G v (t — s,y — z)dz 



+ — \x - y\ U (^-f=4 ) - 1 ) + 2A V t G 2v (t, x-y), (2.20) 



'2vt 



for all t > 0, x,y G E, where f(s,z) = \\u(s,z)\\2 is defined in (2.19). 

Corollary 2.3 (Lebesgue initial data). Suppose that |p(w)| 2 = A 2 (^ 2 +u 2 ) and \i is Lebesgue 
measure. Then for all t > and i,i/Gl, 



E [u(t, x)u (t, y)]=l + (l + O exp 



A 4 t-2A 2 |a;-y| 
Av 

xertc( lX - yl - X2t )-ertc(^) |. (2.21) 



In particular, when y = x, we have 

E[Kt,x)| 2 ] =l + (l + ^ 2 )H(t). (2.22) 

Remark 2.4. If p(w) = u (i.e., A = 1 and <j = 0), then the second moment formula (2.22) 
recovers, in the case n = 2, the moment formulae of Bertini and Cancrini [1, Theorem 2.6]: 



r, / (n(n 2 -l) ) f n(n 2 -l) 

E [\u(t, x) H = 2 exp | | $ ^ _A_J t 

As for the two-point correlation function, Bertini and Cancrini [1, Corollary 2.4] gave the 
following integral form: 



This integral can be evaluated explicitly and equals 

ft — 2\x — y\\ (\x — y\—t 
exp erfc 



Av / V V4z4 



so their formula differs from (2.21). The difference is a term erf (\x — y\/ \jAut). By letting 
x = y in the two-point correlation function, both results do give the correct second moment 
(the difference term is zero for x = y). However, for x ^ y, this is not the case. For instance, 
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as t tends to zero, the correlation function should have a limit equal to one, while (2.23) has 
limit zero. The argument in [1] should be modified as follows (we use the notations in their 
paper): (4.6) on p. 1398 should be 



exp 



/2u 



P € (da)ES 



exp 



The extra term is the last term, which is 



P(T C > t) 



lei 



V2 



exp 



TCS" 



2s 



ds = erf 



V2D 



+ PiTf. > t) . 



erf 



\x — x 



With this term, (2.21) is recovered. 



Example 2.5 (Higher moments for Lebesgue initial data). Suppose that fi(dx) = dx. 
Clearly, J (t,x) = 1. By the above bound (2.15), we have 



E[\u(t, x) | p ] < 2 p - 1 + 2 P/2 - 1 (2 + f) p/2 exp 



8u 



$ I a z - L 2 z 2 




p/2 



We can replace z p by and a Pi? by 2. Then the upper Lyapunov exponent of order p 

defined in (1.5) is bounded by 



X P < 



2 5 p 3 L A 



v 



If q — 0, we can replace a Pi? by \/2 instead of 2, which gives a slightly better bound \ p < 

4 
P - 

A p < 2 3 p 3 X i /u 



2 3 p 3 iA j v. In particular, for the parabolic Anderson model p(u) = \u, we have 



which is consistent with Bertini and Cancrini's formulae \ p = ^p(p 2 — 1) (see [1, (2.40)]). 

Corollary 2.6 (Dirac delta initial data). Suppose that \p{u)\ 2 = X 2 (q 2 +u 2 ) and fi is the 
Dirac delta measure with a unit mass at zero. Then for all t > and x, y G R ; 



E [u{t, x)u (t, y)] = G v {t, x)G v (t, y) - q 2 erfc 

x exp 

In addition, when y = x, we have 

1 



\x - y\ 



2\/vt 
\H-2\ 2 \x-y 
4^ 



+ (— G v/2 ( t 



x + y 



erfc 



x — y\ — \ 2 t 



2Vut 



E[\u(t,x)\ 2 ] = —IC(t,x)+q 2 n(t) . 



. (2.24) 



(2.25) 
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Remark 2.7. If p(u) = u (i.e., A = 1 and <j = 0), then the second moment formula (2.25) 
recovers the result by Bertini and Cancrini [1, (2.27)]: 



E[K*,x)| 2 ] 



2-KUt 



-e vt 




, 7T t t 
1 + \l — 




2v 



which equals /C (t, x; v/2, 1/V Airv) . As for the two-point correlation function, Bertini and 
Cancrini [1, Corollary 2.5] gave the following integral form: 



E [u(t, x)u (t, ; 



2nvt 



exp 



x 2 + y 2 
2vt 



\x - y\ 



1 



o 



y/A-KVt a/s 3 (1 - S 

, , /7rt(l - s) f t 1 - s 
This integral can be evaluated explicitly, and is equal to 



exp 



(x — y) 2 1 — s 
Aut s 



2v 



ds . (2.26) 



G„(t,x)G v {t,y) + ^G* U^ 1 



exp 



i — 2|x — y| 
4^ 



erfc 



|x — y| — t 
\/Aui 



This coincides with our result (2.24) for q = and A = 1. 

Example 2.8 (Higher moments for delta initial data). Suppose that p = 5 and ? = 0. Let 

p > 2 be an even integer. Clearly, J (t,x) = G„(t,x). Then, by (2.15), we have that 

p/2 



E [|u(i, x) | p ] < 2 p - 1 GP(t, x) + 2^' 2 (2G 2 V * £ p ) (t, x 
< 2 p - l G p u (t, x) + 2^- 2 )l 2 L~ p z p p ]C p (t, x 

f y/Anvt 



p/2 



2 p - 1 G p (t,x)+2 p - 1 G p J / 2 2 (t,x) 



I + z 2 L 2 „ 4^ 



7/ 



P P 



P/2 



(the second inequality requires a proof). Hence, for all i6K, the upper Lyapunov exponent 
(1.5) of order p is bounded by 

^ < L p Z P P < z p L p 

p ~ 2v ~ v 

where the last inequality is due to the fact that z p < 2^/p for all p > 2. Note that this upper 
bound is identical to the case of Lebesgue initial data. We can also calculate the exponential 
growth indices explicitly in this case: 



lim - sup logE [\u(t, x)\ p ] < 



a 2 p , L 4 p z 4 



t->+oo t 



\x\>at 



2v 



+ 



2v 



for all a > . 
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Hence, the upper growth indices of order p is bounded by X(p) < zi I?. Similarly, one can 
derive that A(2) > l 2 p /2. Finally, since A(2) < X(p) for all p > 2, we have that, for all even 
integers p > 2, 

\ < Hp) < Hp) < 4 % ■ 

Similar bounds are obtained for more general initial data: see Theorem 2.10 below. 

This following proposition shows that initial data cannot be extended beyond measures. 

Proposition 2.9. Suppose that the initial data is /i = 5' , the derivative of the Dirac delta 
measure at zero. Then the parabolic Anderson model p(u) = Xu (X ^ 0) does not have a 
random field solution in the sense of Definition 2.1. 

2.3 Exponential Growth Indices 

As an application of the above second moment formula, we partially answer the first open 
problem proposed by Conus and Khoshnevisan in [4]: the limits over t in the definitions of 
these two indices do exist when n = 2 and the lower and upper growth indices of order 2 
(see (1.7) and (1.8)) coincide. 

Before stating the main result, we first give some explanation concerning the exponential 
growth indices defined in (1.7) and (1.8). When the initial data is localized, for example, 
when it has compact support, we expect that the position of high peaks of the solution will 
exhibit a certain wave propagation phenomenon. As shown in Figure 1, when a is sufficiently 
large, it is likely that there is no high peaks outside of the space-time cone — the shaded 
region. Hence, the limit over t should be negative. The largest a such that this limit remains 
negative is then defined to be the upper growth index X(p). On the other hand, when a is 
very small, say a = 0, then there must be some high peaks in the shaded region so that the 
limit becomes positive. Hence, the smallest a such that this limit is positive is defined to be 
the lower growth index X(p). 



t t 




Figure 1: Illustration of the exponential growth indices. The initial data, depicted by the 
curve, is localized around the origin. 

Theorem 2.10 (Exponential growth indices). The following bounds hold: 



13 



(1) If \p(u)\ 2 < L 2 (<j 2 +u 2 ) with q = (which implies <£ = q = 0) and the initial data 
H G Mq(M) for some (3 > 0, then for all p>2, 



A(p)< 



a -A 7~ 4 ? 2 T 2 

P ^ + -f±, if o<P< z -^, 

2 2v a v 

z 2 T 2 

7 2 T 2 if R > lph 9 



where z m , m G N, m > 2, are the universal constants in the Burkholder- Davis- Gundy 
inequality. In particular, for p = 2, 

A(2) < \ 2 a 2 " (2.27) 

(2J // |p(w)| 2 > (s 2 +u 2 ) with q = 0, then 

I 2 

A(p) >-£, f or al1 V e -M+(M) ; fi^O and allp>2; 
otherwise, if g ^ 0, i/ien 

A(p) = A(p) = +oo, /or a// /i G .M + (R) and p > 2 ; 

(^) In particular, for the quasi-linear case \p{u)\ 2 = X 2 (q 2 +u 2 ) with A ^ 0, if q = and 
P > #»en 

A(2) = A(2) = A 2 /2, /or a// /i G M^ + (R), fi ^ ; 
otherwise, if q ^ 0, i/ien 

A(p) = A(p) = +oo, /or a// G .M + (R) and p > 2 . 

This theorem generalizes the results by Conus and Khoshnevisan [4] in several aspects: 
(i) more general initial data are allowed; (ii) both non trivial upper bound and lower bounds 
are given (compare with Theorem 1.1 [4]) for the Laplace operator case; (iii) for the parabolic 
Anderson model, the exact transition is proved (see Theorem 1.3 and the first open problem 
in [4]) for n = 2 and the Laplace operator case; (iv) our discussions above cover the case 
p(0) ^ 0. 



14 



Example 2.11 (Dirac delta initial data). Suppose that q = £ = 0. Clearly, 5$ & M. 13 
for all B > 0. Hence, the above theorem implies that for all even integers k > 2, 



G. 



< A(Jfe) < A(fc) < zj^ 2 



This recovers the previous calculation in Example 2.8. 

Proposition 2.12. Consider the parabolic Anderson model p(u) = Aw, A ^ with the initial 
data p(dx) = e~^' x 'dx (B > 0). Then we have 



A(2) = A(2) = < 



— + — ^— «/ < 3 < — 
2 8(3v J ~ 2v 



A^ 
I 2 



2z/ 



This proposition shows that for all 3 G ]0,+oo], the exact phase transition occurs, and 
hence our upper bounds (2.27) in Theorem 2.10 for the upper growth index A(2) are sharp. 



2.4 Sample Path Regularity 

Theorem 2.13. Suppose that p is Lipschitz continuous. Then the solution u(t, x) = Jo(t, x)+ 
I(t,x) to (1.1) has the following sample path regularity: 

(1) If the initial data p is an a-Holder continuous function (a G ]0, 1]) over R satisfying 
(1.4), then 



J G Ci Q (R + xl) U Ci A (R* + x R) , and I G Ci_ i_ (I 
T/iere/ore, 



x 



a.s. 



m = J + I G Ci_ 



X 



U Ci i_ (R+ x R) , a.s. 



f^j If the initial data p is a continuous function satisfying (1.4), iaen 
J G Ci i (R* xl) , and I e Ci_ i_ (R+ x R) , 

Therefore, 



a.s. 



It = Jn + I G C. 



i_i_(RIxR), a.s. 

4 , 2 V t" / 
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(3) If the initial data /i is a signed Borel measure satisfying (1.4), then 
J G Ci A (R* + X R) , and I G Ci_ i_ (R* x R) , 

27iere/bre ; 



a.s. 



u = J + 1 e Ci_ i_ fRI x R) 

4 , 2 \ 



a.s. 



Example 2.14 (Dirac delta initial data). Suppose p(u) = Xu with A ^ 0. If \i = 5o, then 
neither J Q (0,x) nor lim i _ s . 0+ 11^(^^)112 i s continuous in x. For J (0,x) = 5o(x), this is clear. 
As for lim^c^ ||J(t,x)|| 2 , by Corollary 2.6 (with q = 0), we have 

\\I(t,x)\\l = ±1C(t,x)-Gl(t,x) = ^e^$ (^ 2 \f]^j G„ /2 (t,x) . 

Therefore, 

t nm||/( t ,,)|||=f» I'* I- 

<^0+ - +oO if X — . 
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